Quantum-critical behavior of itinerant ferromagnets 
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We study the stabihty of the Quantum Critical Point (QCP) for itinerant ferromagnets commonly 
described by the Hertz-Millis-Moriya (HMM) theory. We argue that in D < 3, long-range spatial 
correlations associated with the Landau damping of the order parameter field generate a universal 
negative, non-analytic Igl^^"''^''''^ contribution to the static magnetic susceptibility ^^(^lO), which 
makes HMM theory unstable. We argue that the actual transition is either towards incommensurate 
ordering, or first order. We also show that singular corrections are specific to the spin problem, while 
charge susceptibility remains analytic at criticality. 
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The critical behavior of itinerant paramagnets at the 
onset of ferromagnetic ordering is the subject of intensive 
experimental (Ij and theoretical 0,13 studies. Of partic- 
ular interest is the existence of a ferromagnetic second 
order transition at low temperatures, and the emergence 
of superconductivity near the phase boundary. The crit- 
ical theory of itinerant ferromagnets was developed by 
Hertz, Millis, and Moriya (HMM) and is based on the 
action of the form: 



S 



d'^qduj„ 



(1) 

where </> is a bosonic field associated with the order para- 
mater, and f is the correlation length, which diverges 
at the phase transition. The two assumptions behind 
Eq. (^) are that the 62n,n>2 terms are non-singular and 
can be approximated by constants, and that the static 
spin susceptibility has a regular momentum depen- 
dence. The HMM theory has been successfully applied 
to explain quantum-critical behavior in a number of ma- 
terials 0, however its key assumptions have been ques- 
tioned in recent studies 0, . The analyticity of the 62™ 
terms was analyzed in detail for an anti-ferromagnetic 
transition where it was demonstrated that, for _D < 3, all 
&2n pre-factors do have non-analytic pieces which depend 
on the ratio between momenta and frequencies of the 
fields Q- For a 2D anti- ferromagnetic quantum critical 
theory (2D fermions and 2D spin fluctuations) these non- 
analytic terms give rise to singular vertex corrections Q . 
For a 2D ferromagnetic case, non-analytic terms in &2n 
are still present, however we found that they do not give 
rise to anomalous exponent in the spin susceptibility for 
D > 1 and therefore are not dangerous. 

In this communication, we question another key as- 
sumption of the HMM theory, namely that the momen- 
tum dependence of the static spin propagator is analytic 
at small q. This assumption is based on the belief that in 
itinerant ferromagnets the q dependence of the (jp term 
comes solely from fermions with high energies, of order 
of Ep, in which case the expansion in powers of {q/ppY 



should generally hold for Pf- This reasoning was 
disputed in Refs. ^ and 6]. These authors considered 
a static spin susceptibility Xsil) in £^ weakly interacting 
Fermi- liquid, i.e., well away from a quantum ferromag- 
netic transition, and argued that for D < 3 and arbitrary 
small interaction, the small q expansion of Xa (<?) begins 
with a non-analytic I^P^^ term. This non-analyticity 
originates from a 2'pp singularity in the particle-hole po- 
larization bubble [a, 0j [3 ^^'^ comes from low- energy 
fermions with energies of the order of vpq <ti Ep. In 
ref. 1^, it was further argued that within RPA, the non- 
analyticity in Xs (q) gives rise to the emergence of a non- 
analytic \q\^~^<pqui term in Eq. Furthermore, the 
prefactor of this term turns out to be negative, which 
signals the breakdown of the continuous transition to fer- 
romagnetism. 

The weak point of this reasoning is that within RPA 
one assumes that fermionic excitations remain coherent 
at the QCP. Meanwhile, it is known 0| that approaching 
the QCP, the fermionic effective mass m* diverges as log ^ 
in D = 3 and as ^^--D smaller dimensions. Moreover, 
m/m* appears as a prefactor of the |(7|^^^ term, hence 
the prefactor vanishes at the QCP. Does this imply that 
Eq. 1^ is valid at the transition? Not necessarily, since 
one has to verify explicitly whether or not the divergence 
of m* completely eliminates a non-analyticity in Xs{q)i 
or just makes the non-analytic term weaker than away 
from QCP. If the latter is true, the non-analytic term 
can still be much larger than q^ at small q at criticality. 

In this communication, we report explicit calculations 
which show that the non-analytic term is still present at 
the QCP and accounts for the breakdown of the HMM 
theory. We first consider the problem in the Eliashberg 
approximation and argue that this approximation leads 
to a non- Fermi liquid behavior at the QCP, but the mag- 
netic susceptibility Xsil) remains analytic in q. We then 
demonstrate that the corrections to the Eliashberg the- 
ory are singular and give rise a universal contribution to 
X7^{(l) in the form q^, with /3 < 2 for Z) < 3 and a 
negative pre-factor. 



The starting point for our theory is a low-energy ef- 
fective spin-fermion Hamiltonian obtained by integrating 
the fermions with energies between fermionic bandwidth 
W and WA < W out of the partition function 
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H = ^Ui.(fc-fcj.)4_^Cfc,„-|-^Xs,d('7)SgS_ 
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FIG. 1: The two lowest-order diagrams for the corrections to 
the static spin susceptibility. The self-energy insertion is E(p), 
the frequency dependent piece 'E{ujm) is already incorporated 
into "zero-order" Eliashberg theory 



where Sq describe collective bosonic degrees of freedom in 
the spin channel, g is the residual spin-fermion coupling, 
and A (the upper cutoff) is roughly the scale up to which 
fermionic dispersion can be linearized near the Fermi sur- 
face. The bare spin susceptibility Xsfiil) is assumed to 
be analytic at g <C p_f and have an Ornstein-Zernike form 
Xs.o = Xo/(?~^ + '?^)- The coupling and xo appear in the 
theory only in combination g = g^xo- 

The perturbation theory for Eq. ^ is an expansion 
into two parameters 3 : 



and X='-^-l = ^ — |^ 



(3) 

where Ef is the electron's Fermi energy and the pre- 
factor in a is chosen for future convenience. We assume 
that g is small compared to Ef, i.e., a is a small pa- 
rameter 0, 0|. That g <C Ef is in line with the very 
idea of a separation between high-energy and low-energy 
physics. At the same time, near a ferromagnetic QCP, ^ 
diverges, and the dimensionless coupling A (x ^ is large. 
Fortunately, one can solve explicitly the strong cou pling 
problem in A while keeping the zeroth order in a [8|, [iJI . 
This amounts to neglecting the vertex correction Sg/g 
and the momentum dependent piece in the self-energy 
E(p), which are of next order in a. The analogous the- 
ory for electron-phonon interaction is known as Eliash- 
berg theory 0. 

Eliashberg theory The Eliashberg theory for a ferro- 
magnetic quantum critical point has been described in 
detail in the literature d, Q| . The results are presented 
here for D = 2, but their validity holds for general Dleqi. 
The elements of the theory are coupled fermionic and 
bosonic self-energies Tj{ujjn) and Ii{q,bj„i) = '^{uim) 
respectively (G'"^(p, w„) = i(w„i -f - vf{p - 

Pf); Xs{q,(^m) = Xo/{q^ + + n(g,w„))). The self- 
consistent solution yields <i>(wm) = 7|^m|7 where 7 = 
gpF/iTTVp), and 



{pFiY^^mg 



E^ 



(4) 



(327r/3y3x)i/3, and Ef = 

2/3, ,1/3 



where /(O) = 1, f{x > 1) 

PfVfI^. At the QCP, ^ = 00 and Y^iiOm) = (^m 
where ujq = a^Ep is the typical bosonic frequency of our 
problem. Note that ujq <C Ef- The w^^^ dependence of 
S implies that at the QCP, the Fermi liquid description 



is broken down to the lowest energies. At the same time, 
Eliashberg theory reproduces the form of the spin propa- 
gator from Eq. CEJ, Xs(g,w„0 oc (^"^ + g^ + 7|t^m|/g)"^- 
Alternatively speaking, in the Eliashberg approximation, 
the Fermi liquid is destroyed at QCP, but the magnetic 
transition remains continuous, and Xs{q) is analytic. 

Beyond Eliashberg theory The validity of the Eliash- 
berg theory is generally based on Migdal theorem that 
states that vertex correction 5g/g and are small 

when the bosonic velocity is much smaler than vp- 
In our case, a'posterioiri analysis of typical bosonic and 
fermionic momenta which contribute to $(ct;) and 
within the Eliashberg theory shows that for the same 
uj ~ ojo, typical qb ^ (li-^oY^^ ^ ctPF, while typical 
fermionic momenta \k — kp\ ~ luq/vp ~ a^pp are much 
smaller. This implies that for a <C 1, the effective bosonic 
velocity vb ^ otvp is much smaller than vp, i.e., Migdal 
theorem is valid. To verify this, we computed the correc- 
tions to dT,(p)/dp and to the static vertex in the limit of 
zero bosonic momentum and for external fermions at the 
Fermi surface and at zero frequency, and indeed found 
Q,i/2 smallness. 

Still, this does not immediately imply that Xs{q) re- 
mains analytic beyond Eliashberg theory. To verify this, 
we have to go beyond the limit of zero momentum, eval- 
uate the correction to the static susceptibility at a finite 
q and check whether it remains analytic in q. Let's do 
this. Consider e.g., the diagram in Fig. ^) with the in- 
sertion of into the particle- hole bubble. Expanding 
this diagram to order q^ we find that it scales as q^I, 
where 



E'^ 

I CX a — TTTT 



1^1^ J {vpi:,-i\n + ^{n)\f V(l^^l7)'/' 

■ (5) 

where P = ll + Py, S{x) cx J dly/{P + \n\-f/\l\) normal- 
ized to S{0) = 0(1) and S{x > 1) cx x'^^^. The in- 
tegrand for / has a highly degenerate pole in the upper 
half-plane of Ix, at vplx cx |17 + S(fi)|. This degenerate 
pole can be avoided by closing the integration contour in 
the lower half-plane of Ix , in which case the non- vanishing 
part in / comes from the non-analyticity in S{x) which, 
we recall, is the boson propagator integrated along the 
Fermi surface. This physically implies that fermions un- 
dergo forced vibrations at typical bosonic frequencies. 
The issue then is what are the typical bosonic x for the 
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non-analyticity in S{x). By analogy with the electron- 
phonon problem, one could expect that corrections to 
the Eliashberg theory come from the processes in which 
typical bosonic frequencies are near bosonic mass shell. 
Then typical Ix are of order {Q'-fY^^ as in the Eliashberg 
theory, i.e., typical x in S{x) are of order one. Substi- 
tuting this typical x into lO, we obtain / ^ 

i.e. 

small and analytic correction to the static susceptibil- 
ity Sxs{q) oc q^a^'^^. However, a careful examination of 
the scaling function S{x) reveals that it is non-analytic 
already at the smallest x, such that there is another con- 
tribution to / from bosons vibrating at typical Ix com- 
parable to those near fermionic mass shell. Indeed, ex- 
panding S{x) in powers of x at a; ^ 1 we obtain 



S{x) = 1 - 



3\/3 

Stt (7|f7|)2/3 



log[ll 



(6) 



The log Ix term is the most important here - its presence 
implies that S{x) possesses a branch cut along the imag- 
inary axis of Ix down to the smallest Ix- One can easily 
make sure that this logarithmic singularity is the conse- 
quence of the non-analyticity of the Landau damp- 
ing piece in the bosonic propagator, i,.e., of the existence 
of the long-range spatial component of the dynamical sus- 
ceptibility. Substituting the small x form of S{x) into jSJ), 
and evaluating integrals we find a much larger, divergent 
contribution to / which behaves as This implies 

that the \ow-q expansion of dxs is actually non-analytic 
in q and begins as g'^/^. Performing an explicit calcula- 
tion without expanding in q, and combining the results 
from the two diagrams for xil) Fig- fflwe obtain for 
the static spin susceptibility 



Xs{q) = 



Xo 



g2 - 0.17|<7|3/2p^^ 



1/2 



(7) 



Observe that the q'^/^ term in ^ is not small in a. This 
is not surprising, since the only excitations involved are 
those near the fermion mass shell, while a measures how 
soft are the mass shell bosons compared to the mass shell 
fermions. Note in passing that the self-energy insertion 
in the diagram of Fig. ^) is not a double counting, since 
the correction to static susceptibility only comes from 
the momentum-dependent piece in the self-energy. The 
frequency-dependent is already incorporated at the 
"zero-order" level, and it does not affect Xs{q)- 

The non-analyticity of the correction to the suscepti- 
bility and the disappearance of a could be also detected 
by analyzing vertex corrections away from the limit of 
q = and zero fermion frequency u>. Indeed, evaluating 
dg{q,uj)/g at finite q and uj, wc obtain 



X ^ \ 1/2 , / ™3 ^ 

q^ Ef 



(8) 



where -0(2; ^ l,y <^ 1) cx ^Jx <C 1, and '4){x ^ 
l,y < 1) = 1-1- Alogy/y/x, A = 0(1). We see 



that, although ip{x,y) < 1 (i.e. vertex correction is 
small in a^^^), the small q expansion of Sg{q)/g begins 



|q| I log i.e., it is non- 



as a^/2|g|| logoj]/^ 
analytic in q and in w, and the Iglogwj term does not con- 
tain a^/2 This singular term again comes from fermions 
near their own resonance and can be traced back to the 
non-analyticity of the Landau damping term. Substi- 
tuting this Sg{q,uj) into the susceptibility diagram and 
performing computations we find that logw dependence 
makes Sxs{q) nonzero in the static limit (without it, 
the correction would be proportional to ^(lu) and van- 
ish at Lu = 0), and the overall power oi q is \q\ from 
the vertex correction times the ratio of typical co and 
vplx- As typical vpqx ~ ^(i^) oc uj'^^^, this ratio yields 



,1/3 



CX (Ix) 



1/2 



q 



1/2 



hence the overall power in 5xs{q) 



is q^^^, as above. 

What happens with higher-order terms? We analyzed 
higher-order particle-hole insertions into the susceptibil- 
ity bubble and found that they add only small, 0{a^^'^), 
corrections to the g^^^ term in Q since at least one inter- 
nal momentum is near the boson mass shell. Higher-order 
particle-particle (cooperon) insertions give rise to 0(1) 
corrections to the q'^/^ term and in principle should be 
included. However, the non-smallness of particle-particle 
insertions is not specific to our problem and is in fact 
customary in Eliashberg-type theories 0, 0, H S ■ 
Particle-particle vertex corrections are known to give rise 
to a pairing instability close to the ferromagnetic QCP 
at Tc = 0{ujq), which implies that there exists a dome 
around the QCP where the normal state analysis is in- 
valid. However, a ferromagnetic T^has an order of mag- 
nitude of roughly Tc < 0.015 loq M Hence the typical 
bosonic momenta for the pairing are of order O.Olapi;', 
much smaller than typical q in Eq. ijJJ). Due to this sepa- 
ration of energy scales, we do not expect particle-particle 
insertions to be relevant. 

Eq. lO is the central result of the paper. It shows 
that the static susceptibility becomes negative around a 
ferromagnetic QCP at g < gi = 0.029 pp- Although 
the pre- factor is small, qi is parametrically larger than 
the typical fermion and boson momenta that contribute 
to the fermion self-energy in the quantum-critical regime 
{\k — kpl ~ ce^PF and qs ~ ctPF, respectively). This 
implies that the whole quantum-critical region is in fact 
unstable. Alternatively speaking, Z = 3 quantum-critical 
behavior at a ferromagnetic transition is internally de- 
stroyed by infrared singularities associated with the non- 
analytic momentum dependence of the Landau damping. 

A negative static susceptibility up to a finite q implies 
either that the instability occurs when is still finite, 
into an incommensurate state with a finite |g| ~ qi |l2l| . 
or the ordering is commensurate, but the magnetic tran- 
sition is of the first order. Which of the two scenarios 
holds remains to be studied. We also verified that away 
from the QCP, when ^ is finite, the q'^^^ term transforms 
into l?!^"^ at the smallest q. The interpolation formula 
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between the (f'l'^ and \q\ forms is rather involved and we 
will not present it here. 

For completeness, we also computed the temperature 
variation of the uniform Xs {q — 0) at the QCP and in 2D 
obtained, with logarithmic accuracy. 





log 


'Ef' 










T 



(9) 



That XsiT,q ~ 0) is negative is another indication that 
the system is unstable at the QCP. 

For arbitrary D < 3, our calculations yield the "correc- 
tion" to in the form \q^^^^^l'^ , again with a negative 
prefactor. In 3D it scales as log \q\. Hence, HMM the- 
ory is destroyed for all D < 3. For D > 3, the correction 
term is smaller than g^, and HMM theory survives. 

Charge susceptibility and gauge theory. We repeated 
the same calculations for the charge susceptibility and 
found in D—2, that the singular |gp/^ terms from the two 
diagrams in Fig. ^ cancel each other, and the behav- 
ior survives. The same holds for the propagator of the 
gauge field - singular Ig^/^ terms from self-energy and 
vertex correction insertions into the particle-hole bub- 
ble again cancel each other. We verified the cancellation 
of non-analytic terms for the second-order diagrams as 
well. Furthermore, we verified in two orders of pertur- 
bation that the momentum dependence of Xc(g, 0) comes 
in powers of {q/ppY, i.e., it is entirely due to the exis- 
tence of a curvature of the dispersion. Physically, the 
distinction between spin and charge susceptibilities is 
in that the ^/\q\ singularity in the dynamical particle- 
hole response function (which generates the q^/"^ non- 
analyticity), is sensitive to a magnetic field, but is in- 
sensitive to the change of the chemical potential. Hence, 
the singularity shows up in the total spin res£onse, but 
does not appear in the total charge response [13 . Math- 
ematically, this distinction is due to the presence, in the 
spin case, of the Pauli matrices in the vertices of the 
diagrams in Fig. such that the self-energy and ver- 
tex correction diagrams do not cancel each other. Note 
that the cancellation or non-cancellation of the singu- 
larities is not directly related to the conservation laws 
for charge and spin susceptibilities, since conservation 
laws only require that Xc(g = 0,Ci;) and Xs(g = 0,Ci;) 
vanish, but impose no constraints on the forms of the 
susceptibilities in the other limit of finite q and zero fre- 
quency. More specifically, Fermi-liquid relations between 
Xs,c(g ^ 0,0; = 0) and XsAl = 0,a; — > 0) imply 
that Xs,c(g — > 0,w = 0) — !■ Cst, but impose no formal 
constraints on the form of the actual q— dependence in 
Xs,c(g,w = 0) at finite q. 

Conclusion To summarize, we studied the stabil- 
ity of the quantum critical point for itinerant ferromag- 
nets, in the limit when the spin-fermion coupling is much 
smaller than the bandwidth, and one would naively ex- 
pect Eliashberg theory to be valid. Within Eliashberg 
theory, the Fermi liquid is destroyed at criticality, but the 



magnetic transition is still continuous, and described by 
bosonic Hertz-Millis-Moriya action. We demonstrated, 
however, that in D < 3, long-range spatial correlations 
associated with the Landau damping of the bosonic order 
parameter field, break the Eliashberg theory and give rise 
to a universal, negative, non-analytic |g|^^+^^/^ contribu- 
tion to the static magnetic susceptibility Xs(giO). This 
term makes the continuous critical theory unstable. We 
argued that the actual transition is either towards incom- 
mensurate ordering, or first order into a commensurate 
state. We also demonstrated that singular corrections are 
specific to the spin problem, while charge susceptibility 
remains analytic at criticality. 
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